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Abstract 



We discuss a magnetic black-hole solution to the equations of motion of the 
string-loop-corrected effective action. At the string-tree level, this solution 
is the extremal magnetic black hole described by the "chiral null model." In 
the extremal case, the string-loop correction is constant, and this fact is used 
to analytically solve the loop-corrected equations of motion. In distinction 
to the tree-level solution, the resulting loop-corrected solution has the hori- 
zon at a finite distance from the origin; its location is a function of the loop 
correction. The loop-corrected configuration is compared with a string-tree- 
level non-extremal magnetic black hole solution which also has the horizon at 
a finite distance from the origin. We find that for an appropriate choice of 
free parameters of solutions, the loop-corrected magnetic black hole can be 
approximated by a tree-level non-extremal solution. We compare the thermo- 
dynamic properties of the loop-corrected and non-extremal solutions. 
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I. INTRODUCTION 



At present, string theory is considered the best candidate for a fundamental theory that 
would provide a consistent quantum theory of gravity unified with the other interactions 
|l|J^. In particular, string theory provides a powerful approach to the physics of black 
holes (review and further refs. in 0,3). In this setting a fundamental problem is that of 
understanding how the intrinsically stringy effects modify Einstein gravity. In this paper 
we discuss two of these effects: presence of scalar fields such as the dilaton and the moduli 
fields and higher-genus contributions modifying the tree-level effective action. We focus on 
higher-genus corrections, because string theory, being a theory formulated on a world sheet, 
intrinsically contains string-loop corrections from higher topologies of the world sheet (these 
vanish only for higher supersymmetries > 4), while a' corrections can vanish in certain 
constructions based on conformal field theories and for a large class of backgrounds. 

We consider the one-string-loop (torus topology) corrections to a special class of back- 
grounds: magnetic black holes provided by the "chiral null models" embedded in (com- 
pactified) heterotic string theory f^J^. String-loop effects modify the gauge couplings in the 
effective action. In the general case, these corrections are rather complicated functions of the 
moduli (fields that describe compact dimensions) 0-0. However, for the extremal magnetic 
black hole solution, the loop correction A is constant. Assuming that the loop correction 
A is numerically small, we obtain an analytic solution to the loop-corrected equations of 
motion. The dilaton of the string-tree-level solution increases at small distances. Thus, at 
small distances, the effective gauge coupling e"*^ + A is sensitive to string-loop correction. 
This produces an important physical effect: the horizon of the loop-corrected magnetic black 
hole is shifted from the point r = to a finite distance determined by the loop correction. 

The loop-corrected solution is compared with a non-extremal string-tree-level magnetic 
black hole solution. Since for small loop correction A the effective string coupling e"*^ at 
large distances is much larger than the loop correction, adjusting free parameters of the non- 
extremal solution, the large-r asymptotic of the loop-corrected metric can be made close to 
that of the non-extremal solution. Requiring that the locations of the horizons of both 
solutions are of the same order, it is possible to fix the free parameters of both solutions 
and to approximate string-loop-corrected black hole solution by a tree-level non-extremal 
configuration. 

In section 2 we review the structure of the 4D magnetic black hole in heterotic string 
theory provided by the chiral null model. 

In section 3 we solve the loop-corrected equations of motion. Starting from the extremal 
magnetic black-hole solution at the string-tree level, we obtain an analytic solution to the 
loop-corrected equations of motion. 

In section 4 we find the position of the horizon of the loop-corrected magnetic black 
hole as a function of the loop correction A and the free parameters of solution. We show 
that by adjusting the free parameters of both solutions the loop-corrected solution can be 
approximated by a non-extremal string-tree-level magnetic solution. 

In section 5 we discuss some thermodynamic properties of the loop-corrected magnetic 
black hole. We calculate the Hawking temperature and geometric entropy of the black hole 
and compare them with the corresponding quantities of the non-extremal solution. Finally, 
we comment on possible microscopic derivation of the entropy. 
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II. DYONIC 4D BLACK HOLE IN TOROIDALLY COMPACTIFIED 
HETEROTIC STRING THEORY 

We begin with a brief review of the chiral null models 0J^. The chiral null model 
is a nonlinear 2D cr-model interpreted as a string world-sheet Lagrangian with nontrivial 
backgrounds: 

L = F{x)du [Bv + K{u, x)du + 2Ai{u, x)dx'] + {dj + Bij){x)dx'dx^ + 7^$(a;) . (1) 

An important special case of the chiral null models (|1]) are the chiral null models with curved 
transverse part of the form 



L = F{x)du[ dv + K{x)duj + f{x)k{x)[dx^ + as{x)dx'^][dx'^ + as{x)dx'^] 

+f{x)k'\x)dx'dx' + bs{x){dx^dx' - dx^dx') + 7^$(x), (2) 



7^ ^ \a'^gi^)R^'\ 

Here ) and V = 2t are non-compact space-time coordinates, u = y2 and 

x^ = Hi are compact toroidal coordinates. Written in "4D form", the chiral null model is 

L = (G;, + B'^^){x)dx^dx'' + G^n{x)[dy'^ + Af'^{x)dx^'][dy^ + 

+ A(2(x)(9l/W-9?/"9x^)+7^0(x). (3) 

The 4Z) string-frame background is related to the fields in the Lagrangian (0) as |]TI | 



r^t 4(1)™ 4(1)" D' D ( A\ 

41)" = Al = G^-^Gm^ , = = $ - ^ In A, A = det G^,, (5) 

and their explicit form can be read off from (^: 

A\ = (0, a,), Al = {K-\ 0), B,^ = (0, 6,), B^, = {F, 0). 



Following [T^ we consider a special class of chiral null models where the functions a<j and 6s 



satisfy the equations 



Here / and k~^ are three-dimensional harmonic functions 

dsd'f = 0, dsd'k-^ = 

and the dilaton (j) is 



The four-dimensional Einstein-frame metric is 



g,, = e-2^G;,. (6) 

From the sigma-model perspective, the AD backgrounds of the above functional form are 
solutions to the equations ensuring conformal invariance of the theory. These backgrounds 
have also the space-time interpretation as black holes that can be obtained by solving the 
equations of motion of the 4D low energy effective action which (in the Einstein frame) is 

S = J d'x^g [r - \{d,ct>f - {d,af - {d,^f 



- - {e-'-^'Fl,, ) . (7) 



Here g^j^^ is the (Einstein frame) metric, a and 7 are the moduli related to the radii of the 
T^. The action (0) is a particular case of the general 4D effective action which is obtained 
by dimensional reduction of lOD supergravity action |jll| by compactifying on the product 



of the tori and by truncating to the relevant set of fields. We shall consider the following 
solution to the equations of motion P,p!^: 



ds" = -A{r)dt^ + A-\dr^ + r^d^l), 
A2(r) = FK'^kf-\ 
2(f) = \nFK~^fk~\ 
e^" = FK, 

e^^ = fK (8) 



where F~^, K, f and are harmonic functions. and F(3) correspond to electric fields, 
F(2) and F(4) are magnetic field strengths. 

Following the general discussion 1^-0, the one- loop-corrected AD action has the following 
form: 



S = l d'x^ (i? - 1(9,0)2 - {d.aY - {d,^f 



+ A)i^f3)-^e-2'^(e-* + A)^(4)) (9) 



Here the loop correction A is a function of Gn = fk and G22 = FK and in the general case 
has a rather complicated form. 

Let us consider purely magnetic [F"^ = K = 1) extremal black holes (/^^ = k). In this 
case in (P) cr = 7 = and backgrounds of the chiral null model are expressed by a single 
function (we consider one-center solution) /q: 
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/o(r) = 1 + -, 
r 

A = f = k-' = /o, 
= b^ = P{1 - cos^), 
= ln/o 7 = 0, (10) 

where and are the nonzero components of potentials in spherical coordinates. The 
magnetic field strengths are = F(^3)ij = Sijkd^fo, and in spherical coordinates have 

a single nonzero component F^^ = — Psin?9. For the magnetic extremal solution, the 
components of the internal metric Gu and G22 are constants, and hence perturbatively [] A 
is also a constant which we assume to be numerically small. In this case the effective action 
(§) takes the form 

s = l d'xV^ (i? - ^id,<pf - {d,^f 
- -/-^{e-f + A)F(i) - ie-2^(e-^ + /\)Ff,^ . (11) 



III. SOLUTION OF THE EQUATIONS OF MOTION 



We look for a static spherically-symmetric solution of the field equations of the loop- 
corrected action (|Tl|). The general ansatz for the metric is 



dsl = -e'^df + e^dr^ + e^dQ^ 



(12) 



The field strength Fij is assumed to have the same functional form as in the tree-level case: 
Fij = —Sijkd^fir), where the function /(r) is to be determined from the field equations. In 
spherical coordinates this ansatz yields the only nonzero component 



F-d^ = /'(r)r^ sin ^9. 
The equations of motion resulting from the action ([TT|) are 

DV+^e-ne^-Ff,)+e-^-Ff3)) = 



(13) 



(14) 



(15) 



R 



+ A) 



1 



g^,Fl^^e~''n2[Fl 



(3) 



0. 



(16) 



One- loop corrections are calculated with tree-level expressions. 
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and 



0. 



(17) 



At the tree level 7 = (see (pUD) and F"^-^^ = F^^y In the first order in A, we have 7 = A71. 
Expanding all the expressions in (|1^) to the first order in A, we obtain 



(0) 



71 = 0. 



Here the subscript (0) stands for the tree-level expressions. From this equation it follows that 
7i = is a solution. Multiplying this equation by 71 and integrating over the space-time, 
we obtain 



- J d\^{D^,f = J d'x^(e-^Ffy 



Since -^(1^3) > (see below), this equation can be satisfied only if 71 = 

Thus, in the first order in A, we have 7 = also, and henceforth we set 7 = 0. 
The field strengths have the following nonzero components 



{F'U = (/'(r)r2 sin ^)2e-^ {F^U = U\r)r^fe 



2\2 -fi 



2^2-2^1 



Eq. (p!^ takes the form 



^2 - 2 



0. 



Away from the origin, Eq. (|T5D reduces to 

(v^(e-* + A)/Vi^^^) = 0. 



(18) 



(19) 



(20) 



Since ~ sin-i?, g'^''^ ~ sin"^-;? and F^^ ~ sin^9, Eq. ( pOD is satisfied identically. To 

account for a singularity at the origin (point-like source of magnetic charge), we calculate 
the flux of magnetic field through an arbitrary sphere enclosing the origin 

j ^F = j d(pd-& sin 't}Ftp^ = Const. 

The Einstein equations ( |T6|) (with one index lifted) take the form |13 



+ tC 



+ A)F- 



(21) 



2e- 



-e 
2 



+ A)F^ 



(22) 



e-^(2/i" + 2z/" + /i'' 



fx'X' - v'X' + /xV) + e 



(e-^ + A)F2 = 0. 



(23) 
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Our next aim is to solve the equations of motion to the first order in A. In the leading 
order (A = 0) we have 



^.W = -ln/o, AW = ln/o, /.(°) = In /o + 2 In r, 0(°) = ln/o, F^^^^ = 2g'^ g' 



2(0) _ o„/2 J _ fi_ 

fo' 
(24) 

in the first order in A we look for a solution in the form 

z/ = -ln/o + An, A = ln/o + A/, /i = In /o + 2 Inr + Am, (25) 

= In /o + Alp, 
F^^ = P{l + Ap) sin^. 
Here n, m, I and (j) are unknown functions, p is a number. One also has 

F^ = 2q'\l + Ar), 

where r = 2p — 2m. In Appendix A we present the detailed solution of the equations of 
motion assuming that m = I = —n Although there are more equations than unknown 
functions, it appears that the ansatz is consistent. We obtain the following solution: 

n = -m = -l = + Ao + A/o + J/oln/o (26) 

Jo ^ 

and 

<^ = -A^A + A, + 2p+ (Ai + l)/o + \k In /o. (27) 
Jo Z I 

Here Ai and p are arbitrary constants. The constants in the above solution are partially 
constrained by imposing the conditions that at large r the metric is asymptotic to the 
Minkowski metric and that the asymptotic value of the dilaton is unity: 

A_x + + = 0, 
+ + + ^ + 2p = 0. (28) 



From the relations (ESD it follows that 



A.,=v^\ (29) 



^There are four equations (|T9|) and ([2l|)-(p3|) for four unknown functions m, n, I and 99. Our choice 
corresponds to the requirement that in the first order in A, as in the leading order, /u + = 21nr. 
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IV. HORIZON AND MATCHING OF SOLUTIONS AT SMALL AND LARGE 

DISTANCES 



Let us discuss our solution. The expressions (Pq)-p7|) were obtained by making two ex- 
pansions: (i) in A, assuming that corrections are smaller than the leading-order expressions: 
|ln/o| > |An|, |ln/o| > |A/|, etc., and (ii) by expanding the exponents e^", e^', etc. to 
the first order in A assuming that 1 > \An\, |A/|, . . .. The bounds that define the range 
of validity of our solution result both in constraints on the free parameters Ai and on the 
domain of variation of r. |^ 

Let us consider the Qoq component of the metric. Introducing new variable 



X 



/o(r) 



r + P' 

we rewrite the expression for the metric component Qqq as 



X 



1 + A L4_i(x- 1) 



(30) 



Here we used relations (^) to eliminate Aq. 

The range of validity of the expression (|30|) is determined by the inequalities 



1 > 



A A_,{x-1] 



X 
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2x X 



and 



In - > 

X 



AL4_i(x-l) + Ai(- 



X 



^ ~^ 2x X 



(31) 



(32) 



Sufficient conditions for validity of the the inequalities ( pi| ) and ( |5^ are their validity for 
each separate term in ( pT]) and (|32D. Solving the inequalities, we obtain 



1 > |AA_i| 
X > \AAi\, 



1 > \AAi\ 

X > I A I In 7—7 . 

A 



(33) 



Note that we assume that |A| -C 1. 

Let us extrapolate the expression ( pO] ) to the region x ~ | A| In |^ and look for a zero (a 
would-be horizon) of the function qqq at small x. For small x we obtain the equation 



It is remarkable that solution (|26|)-(27) is obtained in an analytic form for all r for which we can 
use perturbation expansion in A. Usually it is possible to solve such problem only for small and 
large r and to sew the asymptotic solutions. 

^There is another range 1 > Ai > (e — in which there exists a solution of the inequalities. 
However, for Ai in this interval, the range of small x which we are interested in is excluded. 
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This equation has a solution (zero of the time component of the metric 5^00) provided A < 0. 
Introducing new variable y by the relation x = e"^^"^, we transform Eq. (B¥) to the form 



A, where A = 2 (^j^ — A^i^ e^^"^ . For A^ 1 the asymptotic expansion of solution of 



this equation is WA 



, , . In In A ^ /inlnA' 
hi A - \n\nA + - — — + O - — — 
In A \ InA 



An approximate solution to Eq. 



IS 



Ai + ln(^-A_i) / 1 ^ 



(35) 



Let us consider the dilaton. In Eqs. (P^)-(P^) the parameter A enters through the 
combination A + e~'^, where 



X 



l-A(A_,(l-x) + (A, + l)(i-l) + ^lnl 



(36) 



Near the horizon, at ro ^ r, the function e 



is of order | A| In In this region, A + 



and are of the same order. At large r, e"*^ approaches unity and e"*^ ^ A. This 
suggests that we can look for a solution of the field equations in the whole range tq < r < 00, 
neglecting A as compared to , which has a simpler form than the loop-corrected solution 
and at the same time is sufficiently close to the loop-corrected solution, in particular, has 
the horizon shifted from the origin to r > 0. 

Let us consider a non-extremal solution of the field equations (|14D- (HI) at A = 0. The 
metric is ITf 



1 + 



k sinh ji ' 




1--) dr^ + r^dn 
r / 



(37) 



the dilaton and magnetic field strength are 



^It is important that we considered only smooth decreasing functions 500 and e'^. The component 



of the metric gr 



e in the near-horizon region is a large rapidly varying function. 



^In ref. [15|, as well as in related papers cited therein, the action was obtained by dimensional 
reduction from higher dimensions. The 4D vector fields and scalars appear as metric components 
with a 4D and internal indices and with both internal indices respectively. In the case at hand, 
the action (11) contains two scalar fields (j) and 7 and two vector fields. Although one vector field 
comes from the metric and the other from antisymmetric tensor, it can be verified explicitly that 
the action (O) meets all the requirements of [15|. 
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1 + 



k sinh^/i' 



Fii) 



\3) 



k sinh2/i 
2^ 



-£2, 



(38) 



where £2 is the volume of the unit 2-sphere, k is the position of horizon. The extremal 
limit is /c — s> and ^ 00 with fce^'^ fixed. The Qqq component of the metric interpolates 
smoothly between unity at r = 00 and zero at the horizon. Next, we explore a possibility to 
approximate the loop-corrected solution by a more simple tree-level non-extremal magnetic 
black hole. 

To find a non-extremal field configuration which is sufficiently close to the loop-corrected 
one, it is natural to demand that (i) both solutions have equal magnetic charges, (ii) both 
solutions have identical asymptotic behavior, i.e. at large r the leading O(^) terms of the 
metrics and dilatons of both solutions are equal, and (iii) locations of the horizons of both 
solutions are of the same order. Let us discuss these requirements successively. 

Equating the magnetic charges of both solutions, we have 



k 

- sinh 2/i = P(l - |A|p). 



(39) 



The asymptotics of the component Qqo of the loop-corrected metric (other components 
differ from this expression by the sign of the O(-) term) is 



l--[l + |A|(^,-p+l)l + 0(l) 



where we used the relations (|28| ) and ( pQ]) . The asymptotics of the corresponding component 
of the metric of non-extremal solution is 



(2) 

9oo 



1 



P 

r 



1 + 6 + 



+ 0{ 



1 



where 



k sinh^ a 

6 = ^-1. 



Assuming that the non-extremal solution is close to extremality, from Eq. (pQ]) we obtain 



~ 

2P 



\A\p. 



(40) 



Equating the asymptotics of both metrics and substituting (^OD, we have 

k 1 
2PIAI ~ 4' 



A, 



(41) 



Asymptotics of the functions e are 



-0(1) 



-[l-|A|(A+p+|)]+0(- 
r 4 5 



and 



1- -(1 + 6) +0(1) 
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Using (|4y), we have 



and 



e-*"' = l-f[l-|A|(^+P+l)l 



We see that although it is impossible to obtain strict asymptotic equality e""^'^' = e~^^'^\ 
these expressions are very close if ( ^p^^ + p) ^ 1. 

Further specification requires additional assumptions. Assuming that locations of hori- 
zons of both solutions are of the same order, i.e. ~ P| A| In we have 

Ai~ln^, (42) 

and ^ ^ 

sinh^ u ~ -— - In -— - . 

|A| |A| 

For small |A| <ti 1 we have sinh^/i ^ 1, i.e. the non-extremal configuration is close to 
extremality. Taking the free parameter p = 0(1) and noting that In^ ^ 1, we find that 
both asymptotics of the metrics and dilatons are very close to each other. 

Assuming validity of the strict equality, k = ro, and taking p ~ 1, from the equation for 
the horizon (|35|) we obtain 

/i + ln^ 



A; ~ P- 



1 

|A| 



where we neglected the constants of order 0(1) as compared to Substituting for Ai the 
expression (|4lD , we have k = 2P|A| In y^. 



V. THERMODYNAMIC PROPERTIES 

Let us discuss some thermodynamic properties of the loop-corrected configuration. The 
loop-corrected ADM mass calculated from the asymptotics of the the metric is 



M = P 



1 + |A|(-A_i + Ai 



(43) 



To calculate the Hawking temperature we consider the near-horizon region where we take 
r = tq + p"^ and, assuming that tq ^ p"^, expand the metric in powers of p^. The Hawking 
temperature is determined from the requirement that the resulting metric has no conical 
singularity. We obtain 



hh 



47r 1 + 



(44) 



The geometric entropy defined as one fourth of the horizon area is 
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Sbh = irxl + ^) ■ (45) 

Assuming the validity of the second law of thermodynamics, i.e. that the thermodynamic 
relation SEf,h = T6S is approximately satisfied with 6S = Sbh f\ and T = Tbh, we obtain 
SEbh ~ ''"o ~ -P| A| In On the other hand, taking 

6Ebhr-6M = P\A\{-A^i + Ai + ^), (46) 

we obtain that 

(-A_i + Ai) ~ln — 

which is consistent with the estimate (^T]). 

Let us discuss how the thermodynamic properties of the loop-corrected solution match 
those of the tree-level one. The energy, charge and geometric entropy of the non-extremal 
magnetic black hole (|37D -(p8|) are [l^ 



E ~ A;(3 + cosh2/x), 
Q k sinh 2/x, 

Sbh = 7rfc^cosh^/x. (47) 

Assuming that the non-extremal solution is close to extremality, i.e. ^ 1, the excess of 
energy above the extremal limit is 

5E ~ k. (48) 

Supposing that the non-extremal solution is close to the loop-corrected configuration in the 
sense discussed in the preceding section, and substituting in the expression for 5M the 
expressions for Ai and A_i from (P0|), we obtain again that 5M k P\^\ In j^- Finally, 
we can note that the Hawking temperature Tbh = 4^fccoshV geometric entropy 

Sbh = T^k cosh^/i of the non-extremal solution are close to the corresponding loop-corrected 
expressions. 

The entropy of an extremal dyonic black hole with nonzero electric and magnetic charges 
is equal to Sbh = ■^VQiQ2PiP2- In the non-extremal case the same formula for the ge- 
ometric entropy is valid with the charges Pi = k sinh 7^ cosh 7j and Qi = k sinh (5j cosh (5j, 
i = 1, 2, where 7j and 6i are boost parameters [0. Extremal dyonic black hole is obtained 
in the limit /c — > 0, 7^ — 00, 5j ^ cxd, with Pj and Qi held fixed. The solution (|37D-(p8D can 
be considered as a non-extremal counterpart of the extremal magnetic solution of the chiral 
null model with two magnetic charges (p!OD. The entropy of the non-extremal solution can be 
obtained by the same substitution of charges as discussed above. Setting in the expression 



^Note that at the string-tree level, at A = 0, the entropy of magnetic black hole vanishes. 
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for the entropy of the extremal dyonic black hole Qi = Q2 = k and Pi = P2 = k cosh^/i we 
obtain for the entropy of the non-extremal magnetic black hole the following expression : 



Sbh = 7rA;^cosh^/x, 

which is the same as above ( ^7[ ) 

Statistical entropy of the non-extremal 4D dyonic black hole can be calculated using 



the D-brane technique in the near-extremal limit where all the charges are large ||T^. In 
the case at hand, for non-extremal magnetic black hole, we cannot repeat this calculation 
literally. If, however, in the near-extremal limit, — * 0,7 ^ 00 magnetic black hole can be 
substituted with dyonic black hole with the electric charges Qi = Q2 = k and the D-brane 
counting could be extended to this case, it would provide statistical origin of the entropy of 
the non-extremal magnetic black hole. Anyway, the above remarks can be of relevance in 
other approaches to counting of the statistical entropy 



VI. DISCUSSION 

In this paper we have discussed a solution to the equations of motion of the string-loop- 
corrected effective action in the simplest setting: when the moduli fields and therefore the 
string-loop correction to the gauge couplings are constant. As the tree-level solution we 
chose the magnetic black hole given by the chiral null model. 

The main result of our study is that the horizon of the loop-corrected solution is moved 
away from the origin and that the loop-corrected solution can be approximated by a non- 
extremal black hole. It was found that the requirements that both the loop-corrected and 
the non-extremal solutions have the same large-r asymptotics and positions of their horizons 
are of the same order can be fulfilled if the free parameters of both solutions are connected 
in a special way. 

The string-tree-level chiral null model provides a solution to the low energy effective 
action which in a special renormalization scheme receives no a' corrections 0. The solution 
of the one-loop-corrected effective action we considered (§) is no longer expressed in terms 
of harmonic functions and the a'-corrections are present. However, because now the horizon 
is shifted away from the origin, it is possible that the a' corrections are small and can be 
treated perturbatively. 
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APPENDIX A: SOLUTION OF FIELD EQUATIONS 

In the first order in parameter A the Einstein equations ([2T|)-(p3D are 

^" + m'{q' + V l\\<l' + V y^'^ - + \<l"s = -\hq'\ (Al) 
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m' — h H — ) — ip'q — 2 



/ — m 



Q s = -foq , 



(A2) 



m + n + m / — \- n {—q H — ) + ^ q — q s = j^q . 



l2 



(A3) 



Here 

s = l — (f + T = l — ip — 2m + 2p 
We also need the equation for the dilaton (|1^ in the 0(A) order: 

2 1 9 

Lf + -Lf + -{2m + n - I )q + q s = 0. 
r 2 



(A4) 



(A5) 



We look for a solution such that m = I = —n, because in this case, as at the tree level, the 
components of the metric satisfy the relation ga = g^r- Substituting this ansatz in equations 
lT|)-(|A5|), we have 



m +m(-g +-) + (/? - + -g s = --foq , 



(A6) 



and 



II, II /2 £12 

mq +(fq -q s = foq , 



^" + + q'^s = 0, 
r 



(A7) 



(A8) 



where 

s = 2p — m — if 

Eqs. ( |A2| ) and ( [A3| ) reduce to the same Eq. (|A7|). Combining the equations as ([A6|)- 
(|A7|) + (|A2|) , substituting the expression for s , and introducing new variable u = m + ip, we 
obtain 



u" + u'(g' + -) - q'\u -2p) = 
r 



(A9) 



u' + q'{u- 2p) = foq' 



(AlO) 



and 



p" + -p' -q'^{u-2p) = 



(All) 



Eqs. (|A9|) and ([A10|) are not independent: the first one can be obtained from the second one 

as 



{\m' + -(\m = m 
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Thus, finally we obtain the following system 

u + q'u = ifo + 2p)q' (A12) 

and 

if" + + i2p - u)q'^ = (A13) 
r 

Eq. ( |A12| ) can be integrated yielding 



C , /o , ^ 
Jo ^ 

where C is an arbitrary constant. Using this expression we integrate ( |A13|) : 

^ = C, + CJo + ^ + ^/oln/o (AM) 

Finally, for m = —</? + m we obtain 

m = -C2 + 2p+ (Ci - ^)/o + ^ - ^/oln/o. (A15) 

The relation between the constants used here and those in the main test is: Aq = C^ — 2p; 
= Ci - 1/2; and v4_i = -C/2. 
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